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ABSTRACT 

H^^CN J = 8 — > 7 sub-millimetre line emission produced in the circumstel- 
lar envelope around the extreme carbon star IRC-l-10216 has been imaged at 
sub-arcsecond angular resolution using the SMA. Supplemented by a detailed 
excitation analysis the average fractional abundance of H^^CN in the inner 
wind (^5 X 10^^ cm) is estimated to be about 4 x 10^^, translating into a to- 
tal HCN fractional abundance of 2 x 10~^ using the isotopic ratio ^^C/^'^C = 50. 
Multi-transitional single-dish observations further requires the H^'^CN fractional 
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abundance to remain more or less constant in the envelope out to a radius of 
«i4 X 10^^ cm, where the HCN molecules are effectively destroyed, most proba- 
bly, by photodissociation. The large amount of HCN present in the inner wind 
provides effective line cooling that can dominate over that generated from CO 
line emission. It is also shown that great care needs to be taken in the radia- 
tive transfer modelling where non-local, and non-LTE, effects are important and 
where the radiation field from thermal dust grains plays a major role in exciting 
the HCN molecules. The amount of HCN present in the circumstellar envelope 
around IRC-l-10216 is consistent with predicted photospheric values based on 
equilibrium chemical models and indicates that any non-equilibrium chemistry 
occurring in the extended pulsating atmosphere has no drastic net effect on the 
fractional abundance of HCN molecules that enters the outer envelope. It further 
suggests that few HCN molecules are incorporated into dust grains. 

Subject headings: stars: abundances - stars: AGB and post-AGB - stars: carbon 
~ stars: circumstellar matter - stars: individual (IRC-l-10216) - stars: mass loss 



Introduction 



IRC-l-10216 (CW Leo) has for a long time been used as an archetype for carbon 
stars. The reason for this is the combination of its relative proximity (120 pc) and high 
mass-loss rate (on average 1.5x10"^ yr~^; Schoier & Olofsson 2001) making it possible 
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Aeundez fc Cernichard l2006l ). However, as 



demonstrated by IWoods et al.l (120031 ) in a study of a small sample of high mass-loss rate 
carbon stars, a large spread in molecular fractional abundances, in some cases more than an 
order of magnitude, were found among the sources. 

Surprisingly large amounts of oxygen bearing molecules such as, e. g., H:;.0, H^CO, C^O 



and SiO have been found in the carbon-rich envelope of IRC-l-10216 (IMelnick et al. 
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Haseeawa et al.l l2006l : I Aeundez &: Cernichard l2006l : iTenenbaum et al. 
2006al ). SiO fractional abundance estimates for a large sample of carbon stars also indicates 



that such 'anomalous' chemistries might be a common phenomenon in carbon-rich circum- 
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stellar envelopes (CSEs) (jSchoier et al.l l2006bl ). Various suggestions as to explain the ob 



Fischer- Tropsch catalytic processes, and evaporation of cometary bodies ( 


Bieeine et al. 


2000 


Melnick et al. 


2001; 


Willacy 


2004; 


Cherchnefi 


2006; 


Affundez & Cernicharo 


2006). 



One of the molecules that shows least variation in its fractional abundance between carbon 



stars is HCN, less tha n a factor of two in the sample of IWoods et al.l (120031 ). Furthermore, 
Lindqvist et al.l (120001 ). imaged a handful of carbon stars, with a large spread in mass-loss 
rates, in the HCN J = 1^0 transition using the Plateau de Bure interferometer. Combined 
with a detailed radiative transfer modelling they found that the HCN abundance in their 
sample varied by less than a factor of three in normal carbon stars. A possible explana- 
tion for this is that the HCN abundance is determined close to the photosphere in near 
equilibrium chemistry and that the shocked non-equilibrium chemistry occurring in the ex- 
tended pulsating atmosphere only marginally affects the HCN abundance in a carbon-rich 
environment. 

Nevertheless, HCN is of great importance for the circumstellar chemistry further out 
in the envelope whe re its photodissociation is the start of building mor e complex carbon 



chain molecules (e.g.. lCherchneff et al.lll993l ; ICherchneff &: Glassgoldlll993l ). In addition, line 



emission from HCN has been suggested to be, together with CO, the dominant molecular 
coolant in carbon rich circumstellar envelopes (jCernicharo et al.lll996l ). 



In this paper we present new interferometric observations using the Submillimeter Ar- 
ra}0 (SMA) which, together with multi-transitional single-dish observations, help constrain 
the H^^CN abundance distribution in the circumstellar envelope around IRC-l-10216. The 
observations are described in Sec. 2 and the brightness distribution is discussed in Sect. 3. 
The excitation analysis of the observations is presented in Sect. 4 and the results are discussed 
in Sect. 5. The paper is then concluded in Sect. 6. 



2. Observations 

2.1. The SMA observations 

SMA observations of the H^^qn J = 8^7 line at 690.55207GHz were made on 2004 
April 1 using the compact array configuration with eight 6 m antennas, of which six an- 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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Fig. 1. — Velocity channel maps of H^^CN J = 8 ^ 7 line emission from IRC +10216 obtained 
using the SMA. The contour levels are 15 Jy beam"-"^, where n = —3,-2,2,3,4 (negative 
values have dashed contours), and the beam size is 1'.'6 x O'.'T with a position angle of —58° as 
indicated in the lower right corner of the bottom right panel. The velocity channels (given 
in the LSR frame and indicated in the upper left corner) have been binned to 5 kms~^. The 
systemic velocity is — 26kms~^ (LSR) as determined from CO millimeter line observations. 
Offsets in position are relative to ctaooo = 09''47'"57!39, ^2000 = 13°16'43'.'9. 
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tennas were equ ipped with rece ivers for 690 GHz band observations. Details of the SMA 



are described by IHo et al.l (120041 ) . The adopted source position for IRC+10216 was 0:2000 = 
OgMT'^ST^SQ, ^2000 = 13°16'43'.'9. The weather was good, with a 225 GHz atmospheric opac- 
ity of 0.05 (measured at the nearby Caltech Submillimeter Observatory) and stable through 
the night of observations. The corresponding opacity at 690.55 GHz is 0.98 with a zenith 
transmission of 38%. System temperatures at 690 GHz were in the range from 2000 — 8000 K 
(DSB). The projected baselines ranged from 17 — 68m (39 — 157kA) resulting in a synthe- 
sized beam of 1'.'6 x 0'.'7, with a PA of —58° using natural weighting. In an aperture synthesis 
image the beam size is determined by the diameter of the 'synthesised aperture', which is 
twice the maximum baseline in the case of an uniformly sampled wf-plane. In what follows 
we will adopt the term 'equivalent distance' to describe this transformation from baseline 
separation to angular resolution. Thus, we expect to obtain usable information on scales as 
low as 0'.'7, corresponding to the longest baselines. 

The digital correlator has a bandwidth of 2 GHz and the spectral resolution was 0.406 MHz, 
corresponding to a velocity resolution of 0.18kms~^. Phase and amplitude calibration were 
performed with the Jovian moons Callisto and Ganymede. Bandpass calibration used Cal- 
listo and Ganymede as well as Saturn and Uranus. Observations of Callisto (assuming 120 K) 
provided the flux calibration; the uncertainty of the flux scale is estimated to be ~±30%. 
The data were calibrated using the MIR software package developed originally for the Owens 
Valley Radio Observa tory and adapt ed for the SMA. The calibrated visibility data were im- 
aged and CLEANed (jHogbomI [19741 ) using MIRIAD. Velocity channel maps of the H^^CN 



J = 8 ^ 7 line emission obtained by the SMA are presented in Fig. [H 

The actual analysis and comparison with models are carried out in the wf-plane to 
maximize the sensitivity and resolution of the data. 



Table 1: H^^CN single-dish observations of IRC+10216. 



Transition Frequency Eu Telescope 9^^ rj^ah f T^hdv 

[GHz] [K] ["] [Kkms"'] 

J =1^ 86.340 4 OSO 44 0.60 121.9 

J = 3 ^ 2 259.012 25 JCMT 19 0.69 419.8 

J = 4 ^ 3 345.340 41 JCMT 14 0.63 503.2 

J = 8 ^ 7 690.552 149 CSO 10 0.41 405.0 
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2.2. Single-dish observations 

Multi-transition H^^CN observations of circumstellar emission around lRC+10216 have 
been performed using several telescopes and are summarized in Table [H The individual 
spectra are shown in FigjEl Observations of J = 1 — > line emission were carried out using 
the Onsala 20 m telescopcl (OSO) in February 2007. The Caltech Submillimeter Observatory 
(CSO) was used on March 26, 1994 to detect J = 8^7 line emission. The CSO data are 
particularly important, not only because it probes warmer and more dense material closer to 
the central star, but also since it will provide an estimate on the missing short spacing flux in 
the SMA interferometer data as further discussed in Sect.O In addition, publically available 
J = 3^2 and J = 4 ^ 3 observations using the JCMT telescop^, have been obtained. The 
observed spectra are presented in Sect. 15.11 where direct comparison with model predictions 
are made, and velocity-integrated intensities are reported in Tabled! The data was reduced 
in a standard way, by removing a low order baseline and then binned in order to improve the 
signal-to-noise ratio. The intensity scales are given in main-beam brightness temperature 
scale (Tjnb). 

The OSO and JCMT observations were made in a dual beamswitch mode, where the 
source is alternately placed in the signal and the reference beam, using a beam throw of 
about 11' (OSO) or 2' (JCMT). This method produces very fiat baselines. The CSO ob- 
servations were performed using position switching where off-source data were taken ± 180" 
away in azimuth. The raw spectra are stored in the scale and converted to main-beam 
brightness temperature using Tmb = T^/r/mb- is the antenna temperature corrected for 
atmospheric attenuation using the chopper-wheel method, and r^mb is the main-beam effi- 
ciency. The adopted beam efficiencies, together with the FWHM of the main beam (^mb); for 
all telescopes and frequencies are given in Table [TJ The uncertainty in the absolute intensity 
scale is estimated to be about ±20%. In Tabled] also the energy of the upper level involved 
in the particular transition (Eu) is shown, ranging from 4 K for the J = 1 level up to 149 K 
for the J = 8 level, illustrating the potential of these multi-transition observations to probe 
a large radial range of the CSE (see Sect. 15. ip . 



^The Onsala 20 m telescope is operated by the Swedish National Facility for Radio Astronomy, Onsala 
Space observatory at Chalmers University of technology 

■^Based on observations obtained with the James Clerk Maxwell Telescope, which is operated by the Joint 
Astronomy Centre in Hilo, Hawaii on behalf of the parent organizations PPARC in the United Kingdom, 
the National Research Council of Canada and The Netherlands Organization for Scientific Research 
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3. The H^^CN brightness distribution 

The velocity channel maps of the H^^CN J = 8-^7 line emission obtained by the SMA, 
presented in Fig. [1], show that the brightness distribution towards IRC+10216 appears to 
have an overall circular symmetry. At this level of sensitivity and resolution no signs of 
deviations from a homogeneous wind such as mass loss modulations or clumps are evident. 
There is an offset, on average (O'.'S, -O'.'l), of the H^^QN brightness distribution to the 
adopted source position. The slight offset in the H^^CN map to the nominal source position 
likely results from a combination of residual baseline and passband uncertainties. 

As interferometers lack sensitivity to large scale emission it is of interest to ascertain the 
missing flux of the SMA data. The CSO beam of 10" recovers all the H^^CN J = 8^7 line 
emission from the envelope. Converting the intensity scale in the CSO spectrum from main- 
beam-brightness temperature to Jybeam"^, using a conversion factor of 39JyK~^, about 
840 Jy is recovered by the CSO beam near the systemic velocity (at — 26kms~^ relative to 
the LSR). In comparison, the shortest baselines in the interferometer only pick up about 
Vo of the total emission. 



The CLEANed image of an extended source with missing short-spacing data results in 
an extended depression of negative surface brightness on which the source emission resides. 
This artifact is most evident near the central velocity channels in Fig. [H where the real 
envelope emission is expected to be the most extended. The negative feature will further 
reduce the size and flux of the emission in those channels. Great care needs to be taken 
when analysing the data and comparing to model predictions. For the present set of data 
we find that the most reliable way is to perform the analysis directly in the uf-plane which 
will also maximize the resolution and sensitivity. 

In Fig. [3] the azimuthally-averaged visibility amplitudes are plotted as a function of dis- 
tance to the phase centre in kA. The emission has been averaged over three velocity bins, each 
with a width of 9kms^^. There is a clear trend that the visibilities decrease with increased 
baseline indicating that the emission is well resolved. There is also a significant spread in the 
visibilities, in particular around 50 —100 kA. The innermost region of IRC-l-10216 is k nown to 



be inhomogeneous to some degree (e.g., iMen'shchikov et al.ll200ll : iLeao et al.ll2006l ). While 



it is tempting to assign the deviation of the observed visibilities from a more smooth decline 
with baseline separation to departures from a spherically symmetric and homogeneous wind, 
the quality of the current data set prevents such an interpretation. The implications of the 
spread in visibilities on the modelling, i.e. uncertainties in the derived physical parameters, 
are discussed in Sect. 5.2. 
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Fig. 2. — Multi-transition spectra (histograms) of H^^CN line emission toward IRC+10216. 
Spectra from the best-fit single-dish model (Model 1; solid lines) using a fractional H^^CO 
abundance of 6xl0~^ and a envelope size of 3x10^^ cm are also shown. The hyperfine 
splitting of the J = 1 — > transition significantly broadens the line. This effect is explicitly 
taken into account in the modelling. 
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Fig. 3. — Visibilities, averaged over velocity obtained using SMA interferometer as a function 
of distance to the phase center in kA. (also shown on the upper abscissa is the equivalent 
spatial resolution in arcseconds, see Sect. 12.11 for its definition and meaning). The obser- 
vations are overlayed by model results that have been sampled in the {u, f )-plane in the 
same way as the observations. The best-fit single-dish model (Model 1) using a H^^CN 
fractional abundance /o = 6 x 10"'' and envelope size re = 3 x 10^^ cm is indicated by the 
solid line. Also shown are models using a 50% lower (Model 2; dashed line) and 100% 
higher (Model 3; dash-dotted line) H2 density (corresponding to M = 7.5 xlO^^MQyr"^ and 
M = 3.0 X 10~^ Mq yr~^, respectively) within a radial distance of 1.0 x 10^^ cm. Also shown is 
Model 1 including a temperature decrease of 30% for r < 3 x 10^^ cm (Model 4; dotted line). 
All models are consistent with available single-dish data within Icr. 
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4. Excitation analysis 



4.1. Radiative transfer model 



In order to determine the molecular excitation in the CSE around lRC+10216 a detailed 
non-LTE radiative t ransfer code, based on the Monte Carlo method, was used. The code is 
described in detail in lSchoier &: Olofssoru (120011 ) and has been ben chmarked, to high accuracy , 
against a wide variety of molecular-line radiative transfer codes in Ivan Zadelhoff et al.l (120021 ). 
The CSE around lRC+10216 is assumed to be spherically symmetric, produced by a constant 
mass-loss rate (M), and to expand at a constant velocity (fe). In Sect. 5.2 we will relax the 
assumption of a constant mass-loss rate and allow it to vary in time. 

The physical properties, such as the density, temperature, and kinematic s tructure pre- 

vailing in the circumstellar e n velop e arou nd IRC-l-10216 have be en determined in lSchoier fc Olofsson 
( I2OOOI . I2OOII ) , ISchoier et al.l (120021 ) , and ISchoier et al.l (l2006b[ ) , based on radiative transfer 
modelling of multi-transition CO line observations, from millimetre to IR wavelengths (the 
reader is referred to these works for more specific details of the input parameters). This 
model, where the circumstellar envelope is formed by a mass-loss rate of 1.5 x lO~^M0yr~^ 
and expanding at a velocity of 14.0 km s~^, is used as input to the H^^CN excitation anal- 
ysis. The inner radius of the model is taken to be ro = 1.7 x 10^^ cm (= 11 AU = 0'.'095). 
The kinetic temperature structure is shown in Sect. 15.11 where the molecular excitation is 
discussed. 



4.2. The H^^CN abundance distribution 

The abundance distribution of H^^CN is assumed to be described by a Gaussian 

/W = /o exp (^^y j , (1) 

where / = ?t,(H^^CN)/?t,(H2) is the fractional abundance, i.e., the ratio of the number density 
of H^^CN molecules to that of H2 molecules. Here, /o denotes the photospheric fractional 
abundance of SiO and it is assumed that some process effectively destroys the molecules at 
r > Te such as, e.g., photodissociation by the ambient interstellar uv-field. 

The excitation analysis includes radiative excitation through the stretching mode at 
3 fim and in the bending mode at 14 fim. The stretching mode at 5 fim includes transitions 
that are about 300 times weaker and is therefore not included in the analysis. In each 
of the vibrational levels we include rotational levels up to J = 29. Hyperfine splitting of 
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the rotational levels were included only in the J = 1 levels (where the splitting is larger 
than the local turbulent width) and the resulting line overlaps were accurately treated as 
described in Lindqvist et al. (2000). Also, /-type do ubling in the 1 4 /xm transitions was 
included. Relevant molecular data are summarized in Schoier et al.l (l2005h and are made 
publicly available through the Leiden Atomic and Molecular Database (LAMDA)0. 

For IRC+10216 thermal dust emission provides the main source of infrared photons 
which excite the vibrati onal states. The add ition of a dust component in the Monte Carlo 
scheme is described in ISchoier et al.l ( 12002| ). The dust-temperature structure and dust- 
density profile for IRC-l-10216 ar e obtained from radiat ive transfer modelling of the spectral 
energy dist ribution using Dusty Jllvezic fc Elitzurl 119971 ). Details on the dust modelling can 
be found in lSchoier et al.l (l2006bl ). The best-fit dust model is obtained for an optical depth at 
10 /im of 0.9, a dust-condensation temperature of 1200 K (corresponding to a radial distance 
of ro = 1.7 X 10^^ cm and also taken as the inner radius in the H^^CN excitation analysis) 
and an effective stellar ter nperature of 2000 K. The lumino sity is 9600 L©, obtained from a 
period- luminosity relation (iGroenewegen fc Whitelocklll996l ). and the corresponding distance 
is 120 pc. 



5. Results and discussion 

5.1. The extended envelope 

In order to determine accurate abundances it is important to know the spatial extent 
of the molecules. Models where the size of the emitting region (re) and the photospheric 
abundance (/o) in Eq. [T]are varied simultaneously have been performed. The best fit model 
is found by minimizing the total defined as 

N 

2 _ \^ 

Xtot / ^ 

i=l 

where I is the integrated line intensity and a the uncertainty in the measured value (usually 
dominated by the calibration uncertainty of ±20%), and the summation is done over N 
independent observations (here A^ = 4). 

First the available observed multi-transition single-dish data (see Fig. [2]) are modelled. 
These various rotational transitions require different excitation conditions (see Table [1]) and 



(-/^mod 



obs I 



a 



(2) 
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will probe the CSE at different spatial scales. The energy of the upper level involved in 
the particular transition {Ey) ranges from 4 K for the J = 1 level up to 149 K for the J = 8 
level. This means that the lines will not behave in the same way to changes of the ad- 
justable parameters and /o (see Eq. 1) allowing for the possibility to constrain them. The 
sensitivity of the line emission to variations of and /o is illustrated in Fig. |H We find 
that it is possible to simultaneously determine both the spatial extent of the H^'^CN molec- 
ular envelope, re = 1.0±0.8 x 10^^ cm, and the fractional abundance of H^^CN molecules, 
/o = 7.0 ± 2.5 X 10^''. However, it is hard to put good constraints on the exact H^'^CN enve- 
lope size from modelling of single-dish data alone. 

Interferometric observations can help to further constrain the HCN envelope size. In 
Fig. [3 observation s of H^^CN J = 1 — » lin e emission towards IRC+10216 using the BIMA 



interferometer by iDayal fc BiegingI (119951 ) are presented. The sp atial resolution o f 8^^ is 



enough to constrain the envelope size to within 2.5 — 6.0 x 10^^ cm. iDayal fc BiegingI (119951 ) 
in their excitation analysis found a fractional abundance /o(H^^CN) = 7.8 x 10~^ and an 
envelope size Ve = 2.4 x 10^^ cm, in good agreement with the present analysis after correcting 
for their slightly smaller distance and larger mass loss rate. Lindqvist et al. (2000), in 
their modelling of PdB interferometric HCN J = 1 — > observations of IRC+10216, found an 
envelope size of 4 x 10^^ cm and /o(HCN) = 5 x 10^^ which translates to a fractional H^'^CN 
abundance of 1 x 10~^ using a ^^C/^^C-ratio of 50 (Schoier & Olofsson 2000) in agreement 
with the present analysis. It should be noted that Lindqvist et al. did not properly take 
into account a dust component in their excitation analysis (explaining their somewhat larger 
fractional abundance), instead they introduced a central, cool (510 K), blackbody to simulate 
the dust excitation (see the problems with this simplified strategy discussed by Lindqvist 
et al. 2000). The collisional rate coefficients used were also somewhat diffe rent compared to 



the pr esent analysis. Additional ab undance estimates for IR C+10216 by ICernicharo et al. 



(119961 ) based on ISO/LWS data and ICernicharo et al.l (119991 ) in their analysis of ISO/SWS 
observations found fractional abundances of HCN in the range 1 — 3 x 10~^, again in good 
agreement with the present analysis. 

Model spectra from the best-fit single-dish model, i.e. the model with the lowest value 
in Fig. |l](hereafter Model 1), are directly compared to observations in Fig. [21 The quality 
of the best-fit model, where /o = 6xlO~'' and re = 3xl0^^cm, is exceptionally good with a 
reduced (= Xtot/ (^ " 2)) of 0.1. Also the extra line broadening in the J = 1^0 transition, 
due to the hyperfine splitting, is very well reproduced. Moreover, the increased amount 
of self-absorption due to the higher optical depths with increasing rotational transitions, 
producing more triangular shaped line profiles and moving the peak of the line towards the 
red part of the spectrum, is well reproduced by the model. 
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Fig. 4. — x^-map showing the quahty of the fit to available single-dish data when varying 
the adjustable parameters, /o and rg, in the model. Contours are drawn at the 1, 2, and 3 a 
levels. Indicated in the upper right corner are the reduced of fhs best-fit model and the 
number of observational constrains used, N . 
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Fig. 5. — Azimuthally averaged (near the systemic velocity) H^^CN J = 1^0 interfer- 
ometric BIMA observations (Dayal & Bieging 1995) toward IRC+10216. Models using 
/o(Hi3CN)=9xlO-^ and re = 2.5x10^^ cm(dash-dotted line) and /o(Hi^CN) = 7x10"^ and 
Te = 6.0x10^^ cm (dashed line) are also shown. The FWHM (8") of the clean beam is indi- 
cated. 
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Results from the excitation analysis are shown in Fig. [6] for Model 1 where the excitation 
temperature, Tgx, and the tangential optical depth, rtan are plotted as functions of radial 
distance from the star for the J =1-^0 and J = 8 — > 7 transitions. It is clear that both lines 
are formed in non-LTE conditions far from the kinetic temperature of the gas (Tkin, 

also shown in Fig. [6]) where rtan is peaking. The J =1^0 line emission (shown here for the 
inherently strongest of the three hyperfine components the F = 2^1 transition) is optically 
thin where as the J = 8 — > 7 emission is moderately thick. From the excitation analysis it is 
also clear that the J =1^0 line emission is mainly produced in the external cooler layers 
of the H^'^CN envelope and it is thus more sensitive to changes in Ve than the J = 8^7 
transition which is more sensitive to /q. As an example, reducing the envelope size in 
Model 1 by 50% the J = 1 ^ line emission, as observed by the OSO beam, gets 49% weaker 
whereas the J = 8 ^ 7 line emission, in the CSO beam, is lowered by only 4%. If instead the 
fractional abundance in Model 1 is lowered by 50%, line emission from the J=l—>-0 and 
J = 8^7 transitions are reduced by 34% and 31%, respectively. In all, this illustrates the 
usefulness of multi-transition observations in lifting the degeneracy between envelope size 
and fractional abundance which otherwise severely hampers an excitation analysis when no 
direct information of the spatial distribution of the molecules exist. 

In the excitation analysis it is important to include also vibrationally excited states. 
Retaining only the ground vibrational state in Model 1 lowers the line intensities of the 
observed transitions by 42% on average. This means that in a multi-transitional analysis, 
since the various transitions behave in a different manner, not only the derived fractional 
abundance will be severely affected but also, to a lesser extent, the physical size of the H^'^CN 
envelope. 

5.2. The inner wind 

The interferometric SMA observations of H^'^CN J = 8^7 line emission provide addi- 
tional constraints on the H^^CN modelling. From the excitation properties of the J = 8—^7 
transition discussed in Sect. 15.11 it is clear that this line will mainly provide information 
on the warmer more dense circumstellar material closer to the star, typically within 5" in 
radius. This is further accentuated by the lack of short spacing information, and hence loss 
of sensitivity to extended emission, in the set of data analysed here (see Sect. 12. ip . 

The visibilities as a function of baseline shown in Fig. [3] should be regarded as the best 
constraint on the model as it is then tested over a large range of spatial scales. The inversion 
process from the [u, f )-plane to the image domain will also suffer from artefacts introduced 
by the large amount of missing flux, and although this should be the case also when analysing 
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Fig. 6. — Results for the best-fit single-dish model (Model 1) from the excitation analysis, 
a) The kinetic temperature of the gas particles as derived from the CO modelling is shown 
by the solid line. The dotted and dashed lines gives the excitation temperature of the H^^CN 
J=l—>-0, F = 2^1 and J = 8^7 transitions, respectively, b) The dotted and dashed lines 
gives the tangential optical depth (rtan) at the line center of the H^^CN J = 1^0, F = 2— >1 
and J = 8— s>7 transitions, respectively. 
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the model results, this may enhance deviations from the model. 

Based on the best-fit model from the single-dish analysis a brightness map of the emis- 
sion around IRC+10216 was produced. To this map the same spatial filtering as present 
in the observations was applied and azimuthal averages produced in the {u, f )-plane. The 
observed visibilities can then be directly compared with those from the model. As shown 
in Fig. [3] the best-fit single-dish model (Model 1; dotted line) provides a reasonable fit to 
the observed visibilities over the full range of baseline separations, both for velocities near 
the center and close to the edges of the line. This also means that the data is consistent 
with a constant expansion at the terminal velocity (14kms~^) of the wind. The spread in 
visibilities that is observed can be accounted for by modulating the mass-loss rate. In Fig. [3] 
models with a 50% lower (Model 2; dashed line) and 100% higher (Model 3; dash-dotted 
line) H2 density (corresponding to M = 7.5 xlO~^M0yr~^ and M = 3.0 xlO~^M0yr~^, re- 
spectively), within a radial distance of 1.0 x 10^^ cm are shown. The fits to the single-dish 
data are still within the la limit. With the uncertainties of the observations in mind, we 
find the single-dish data and the SMA interferometric data to be consistent with an H^^CN 
abundance distribution described by Eq. [1] with /o = 4 x 10~^ and re = 4 x 10^^ cm. 

The SMA data can not provide an estimate of the H^^CN envelope size, within reason- 
able limits, due to the notorious problem of filtering out of large scale weak emission that 
hampers interferometers and the fact that the J = 8^7 emission is mostly excitation limited 
rather than photodissociation limited (see Sect. 5.1). 



5.3. HCN line cooling 



It has been suggested that HCN line emission coul d be an important coolant in the inner 
parts of circumstellar envelopes around carbon stars (iGroenewegenl 1 19941 : ICernicharo et al. 
19961 ). In Fig. 15.21 we show the estimated total HCN line cooling from our best-fit single dish 
model (Model 1; solid line). The line cooling from H^^CN can be self-consistently calculated 
once the excitation is known from 



(3) 



j,>i 



where Ui denotes the level populations of H "^CN, •jui the coUisional rate coefficients, rin^ is 
the H2 number density, and C is the emissivity per molecule in ergs~^cm^^. In order to 
obtain the total HCN cooling this result was multiplied by the isotopic ratio ^^C/^'^C = 50. 
Within a radial distance of about 2 x 10^^ cm 1") HCN line emission dominates over 
that from CO. Thus a lower kinetic temperature is expected. To simulate this effect the 



- 16 - 



kinetic temperature structure was lowered by 30% in this region. This has no effect for the 
single-dish data given their larger beam-sizes (> 10") but lowers the intensity of the H^'^CN 
J = 8—>-7 line emission from the region where HCN line cooling is effective, still providing a 
reasonable fit to the observations (Fig. [31 dotted line; Model 4). 



5.4. Comparison with chemical models 



Recent s t udies o f SiO line emission around lRC+10216 by ISchoier et al.l (l2006al ) and 
Schoier et al.l (j2006bl ) have shown that non-equilibrium chemical processes are required in 
order to explain the very high amounts (almost two orders of magnitudes higher than ex- 
pected from an equilibrium chemistry) of SiO present in or close to the photosphere. In 
addition, at ~ 5 stellar radii there is a drastic drop in the fractional abundance by about an 
order of magnitude explained by effective adsorption onto dust grains. The effect of shocks 
is strongly dependent on the C/O-ratio and in a carbon star such as lRC+10216 the minor 
species, such as SiO, are expected to be much more strongly affected than the most abun- 
dant ones, such as HCN. Also, SiO is very refractory compared to HCN, which means that 
SiO will condense at higher temperatures much closer to the star, where the densities are 
significantly higher, with the result that the fractional abundance of HCN in the gas phase 
is expected to be much less affected by this process. 

LTE stellar atmosphere models pr edict that the HCN fractional abundance in carbon 



stars is in the range 1 — 5 x 10 ^ (e.g.. lWillacy fc Cherchnefflll998l : ICherchneffll2006l ) which 



agrees well with the value of approximately 2 x 10 ^ derived for IRC +10216 in the present 



analysis and with circumstellar values for carbon stars in general (ILindqvist et al.l 12000 



Woods et al.l 120031 ). This would suggest that HCN is formed in the photosphere under near 
LTE conditions. 

Recent non-equilibrium chemistry models of the dynamically pulsating atmosphere of 
AGB stars indicate that the passage of shock waves can have a significant effect on the 
abundance of HCN close to the photosphere. It is likely that both the C/O-ratio and the 
shock parameters can have a large effect on the processing of HCN in the innermost part of the 
wind (1 — 5 stellar radii) . The fractional abundance of HCN just before dust formation sets 
in is estimated to be f» 4 — 6 x 10^^ for carbon stars (jWillacy &: Cherchnefflll998l : ICherchneff 
20061 ). i.e., up to an order of magnitude lower than those predicted from equilibrium models. 
For an extreme carbon star such as IRC-l-10216, where the C/O-ratio is ~1.5, the initially 
high HCN fractional abundance (relative to H2) gets reduced to about 4.4 x 10~^ at 5 stellar 
radii (before dust condensation), i.e., a factor of about five lower than the circumstellar 
value obtained in the present analysis. However, given the uncertainties involved in both 
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the chemical models and the derived HCN fractional abundance from observations we find 
that it is not possible to rule out the importance of non-equilibrium chemical processes in 
the formation of HCN in the (extended) atmospheres of IRC+10216. Possibly, a statistical 
study based on a large survey of circumstellar HCN line emission, from AGB stars with 
varying photospheric (e.g., C/O-ratio) and circumstellar (e.g., mass- loss rate) properties, 
could provide better constraints on the chemical models. Such work is currently underway 
(Schoier et al., in prep.). The apparent consistency between LTE stellar atmosphere models 
and the derived circumstellar HCN fractional abundance in IRC-l-10216, together with the 
suggestion that non-equilibrium chemical processes in the extended atmosphere will serve 
to reduce the photospheric abundance, indicates that few HCN molecules are incorporated 
into dust grains. 

6. Conclusions 

A detailed excitation analysis, based around a non-LTE and non-local radiative transfer 
code, has been performed for H^^CN in the circumstellar envelope around the extreme car- 
bon star IRC-l-10216. New sub-arcsecond interferometric observations of H^^CN J = 8^7 
line emission, performed using the SMA, put constraints on the HCN abundance distribution 
on scales <5 x 10^^ cm (= 330 AU = 2'.'8). In combination with multi-transition single-dish 
observations, probing the emission on larger spatial scales, it is found that the H^^CN frac- 
tional abundance remains more or less unaffected in the envelope at about 4 x 10"'' until 
photodissociation sets in limiting the size to about 4 x 10^^ cm (=2700AU = 22"). In the 
excitation analysis of HCN it is of utmost importance that radiative excitation through the 
vibrationally excited states is properly accounted for. Also, HCN provides effective hne 
coohng in the inner part of the envelope (< 2 x 10^^ cm) that is higher than that from CO. 
The fractional abundance of HCN derived in the circumstellar envelope of IRC-l-10216 is 
consistent with model atmosphere predictions and suggests that non- equilibrium chemical 
processes in the extended pulsating atmosphere has no drastic net effect in regulating the 
amount of HCN molecules that enters the outer wind. 
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Fig. 7. — HCN line cooling obtained from the best-fit single-dish model (Model 1; solid 
line) compared with that obtained from CO (dashed line). The same underlying physical 
structure has been used in the two cases. 



